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Abstract. Space- VLBI observations of stellar sources represent a challenge 
since there are few sources with sufficiently high brightness temperature for de- 
tection on space-ground baselines. X-ray binaries (XRB) are among the few 
types of stellar radio sources that can be detected. Observations of the unusual 
X-ray and 7-ray binary system LS I +61° 303 obtained with the HALCA satellite 
and a 20-element ground array are described. The data in this 48-hour exper- 
iment represent some of the best quality VLBI observations of LS I +61° 303. 
No fringes were detected on HALCA baselines. 10-minute snapshot images were 
produced from the global ground array data and reveal an expansion velocity 
of 800 km s _1 . Some of these image data reveal hints of more extended emis- 
sion but high-SNR closure phase data do not support relativistic outflow in the 
plane-of-the-sky in LS I +61° 303. The largest closure phase rates are consistent 
with an outflow of ~ 1000 km s _1 as deduced from the image data. The closure 
phases also show no evidence of structure variation on size scales greater than 
~ 10 mas. A number of issues related to VSOP2 observations of stellar radio 
sources are raised. 



1. Introduction 

A challenge of observational astrophysics with radio interferometers is the funda- 
mental limitation the observing array places on brightness temperature sensitiv- 
ity. The sensitivity of the individual telescope elements in the array determines 
the flux density limit. The HALCA spacecraft had a 7a sensitivity in a few 
minutes to a ground-based telescope of ~ 100 mJy, and so for an array with 
B max ~ 10, 000 km (a resolution of ~ 1 mas at 5 GHz) this implies a brightness 
temperature limit of Tb ~ 10 10 K. Consequently the types of stellar sources 
that could be detected on baselines to HALCA was limited to masers, pulsars, 
and some XRBs. There are undoubtedly others types of sources that could be 
detected, but they will exhibit either outburst phenomena or coherent emission 
processes. 

The subjects of maser and pulsar observations are well covered by other pa- 
pers presented at the meeting and the previous VSOP symposium. In this paper, 
the specific case of VSOP observations of the high mass XRB LS I +61° 303 will 
be described. Special attention will be given to some properties of the radio 
emission in LS I +61° 303 as a means to reveal some of the challenges that need 
to be addressed for future space- VLBI observations of stellar sources. 
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What is LS I +61°303? 



Th e luminous early-t y pe st ar LS I +61° 303 was discovered as a radio source 
by IGregorv &; Tavloil (1 19 781). who proposed an association with the COS-B j- 
ray source 2CG 135+01 (|Hermsen et al.l[l977f ). The radio emission was shown to 
have regular, non-thermal outbursts with a period of 26.51 days (jTavlor Gregory 
1984); the flux density rises from a quiescent level of around 20 mJy to a peak of 
up to 300 mJy over ~ 2 days, and then decays over ~ 8 days to its quiescent level. 
The period is related to the highly eccentric orbit of a compact object around 
a rap idly rotating BOe star (jHutchings &: Cramptonl Il98ll ; iGrundstrom et ah 
2007| ). Th e 26.5-day periodic behavior is also manifest at o ther wavelength re- 



gions: IR dParedes et al.lll994h. X -rav dParedes et al.l ["L9971 ) and, most recently 
7-ray energies ([Albert et al.ll2006l ; lMaierll2007f T~ 

Variable TeV emission has been detected in three other high-mass XRB 
systems. Together they form a group of objects that have been referred to as 
Binary TeV sources (BTVs). PSR B1259-63 is a 3.4-yr period massive binary 
syste m with a pulsar in an eccentric orbit around a BOe star ( Johnston et al.l 
1992). The TeV and X-ray emission is interpreted as inverse-Compton (IC) 
emission due to relativistic electrons accelerated in a wi nd-collision between th e 
pulsar wind and that of the massive companion (e.g. Tavani &: Aronsl 1 19971 ). 
Cygnus X-l is a well-established mic r o-quasar system w ith clearly observed rel- 
ativistic jets (e.g. IStirling et aill200ll ; iFender et al.ll2006l ). The situation in both 
LS 5039 and LS I +61° 303 is far from clear. In LS503 9, there is eviden ce of 
relativistic jets but the picture is not completely clear (jRibo et al.l 12008). In 
LS I +61° 303, both micr o-quasar (e.g. iMassi et al.ll2004l ) and pulsar- wind neb- 
ulae ( Dhawan et al.ll2006t) models have be en advanced, with the merits of each 
scenario discussed bv lRomero et all (|2007 , and references therein) . The remain- 
der of this paper will discuss the insights revealed by VSOP observations of 
LS I +61° 303 and how they may inform this discussion. 



3. Previous VLBI observations 

VLBI observations of LS I +61° 303 have been completed by a number of differ- 
ent groups. Observations w ith the EVN near p eak flux indicate a low expansion 
velocity of ~ 600 km s" 1 dTavlor et al.lll992n . However, an expansion veloc- 



ity of ~ 18, 000 km s was measured dur ing a mini-outburst observed during 
the main outburst ( Peracaula et al.l Il998l ). More recently, evidence for rela- 



tivistic expansion has been presented from EVN and MERLIN observations. 
Massi et al.l ( 200ll ) observed an extension of the central source to the SE. The 



lack of a corresponding symmetric feature to the NW led to an interpretation of 
a Doppler amplified outflow close to the line-of-sight with an inferred ejection 
speed of 0.4c. A more recent MERLIN observation suggests a similar structure, 
though with a double-sided jet that extends up to 200 AU ( ~ 70 mas) from th e 
central source and which precesses significantly in 24 hours (jMassi et al.l 120041 ) . 
This is interpreted as due to a jet with an outflow speed of 0.6c. Clearly, these 
observations do not present a consistent picture of the dynamical evolution of 
the radio emission arising in LS I +61° 303. 
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Figure 1. Left: The epoch of VSOP observations (horizontal bar) relative to 
an outburst in LSI +61° 303 as observed by the GBI. The orbital phase from 
ICasares et al.l (|2005l ) is shown, and the phase of periastron. Right: Visibility 
distribution of LS I +61° 303 over 1 hour. The source is clearly resolved. 



4. VSOP observations of LS I +61°303 

Observations over 48 hours were obtained on 1999 September 17 at 5 GHz with 
HALCA and a 20-element ground array including some of the largest telescopes 
available (the phased- VL A, Effelsberg, Jodrell Bank). The goal was to image 
LS I +61° 303 with sub-milliarcsecond resoltion to observe structural changes as 
the source evolved during the 2 days of outburst. Though VSOP observations 
pre-date some of the observations described in Sec. [33 they represent perhaps 
the largest global array assembled to observe LS I +61° 303, and such high- 
SNR data is arguably t he best avail a ble to day for this source. They have been 
described previously bv lTavlor et al.l ( 20001 ) . but more recent analysis is included 



here. 

Unfortunately, the observations started approximately 60 hours after the 
onset of outburst (Fig [I]- left), and data was collected during decline. Further- 
more, fringes on baselines to the HALCA spacecraft were not detected. The 
visibility distribution from part of the observation is shown in Fig [T] (Right). 
This plots demonstrates that at the spatial frequencies sampled by HALCA 
baselines (> 200MA), the correlated flux (~ 20 mJy or less) is too low for direct 
fringe detection (see comments in Sec. |5T|) . 

The ground-telescope array data alone provides a very high-quality VLBI 
observation of LS I +61° 303 during the decline phase of an outburst. Standard 
phase-referencing techniques were used in the ground observations, enabling 
antenna gain calibration using the bright quasar J0228+67 at a few degrees 
distance every 20 minutes. Snapshot images every 10 minute were generated, 
with two examples from 112 total images shown in Fig. [2j 

The image data show the source has a maximum extent of ~ 4 mas, similar 
to previous observations taken during quiescence. Comparison of visibility data 
on each of the two days shows evidence for a ~ 0.3 mas angular expansion 
over 24 hours, correspon ding to an expansion of 0.6 AU at a distance of 2 kpc 
( Frail Sz Hiellminglll99ll ). and an inferred expansion velocity of ~ 800 km s _1 . 
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Figure 2. Two 10-minute images from the 113 in the observation. Contours 
are 8, 12, 18 . . . 80% of peak. Frame 21 (left) and Frame 79 (right) are typical 
examples. Visibility modelling of data from each day reveals expansion of 
~ 0.3 mas over 24 hours. 



A feature seen in some of these images is low surface brightness extensions to 
the central emission region. These extensions are roughly symmetric about the 
central source and have been interpreted as emission from high-velocity outburst 
ejecta. If they are associated with the outburst the inferred expansion velocity 
is 42, 000km s _1 , and broadly consistent with expansion deduced from MERLIN 
and EVN observations by Massi et al. 

However, interpretation of the extensions in some of the images is hindered 
by a) the inability to follow the features from one 10-minute snapshot to the 
next and b) the symmetric nature of the features. This led us to consider if 
these were artefacts of calibration errors. A robust technique for eliminating 
calibration artefacts is to examine closure quantities, particularly phase where 
any time variation in closure phase is due only to intrinsic properties of the radio 
source. This is particularly useful with high-SNR data since even small closure 
phase variations can be readily detected. 

Closure phase for a number of antenna triangles is shown in Fig. [3] and 
several things can be deduced readily. Firstly, on 10-mas size scales represented 
by the EB-JB-MC triangle (bottom frame), the closure phase is identically zero. 
This implies LS I +61° 303 is a point source on these angular scales with no 
variable structure during this observation. The next important feature is that 
the largest changes in closure phase are ~ 20° over a few hours. What does 
this imply? 20° is approximately 1/20 of the beam so for a 1 mas beam this 
corresponds to an asymmetric structure changes of ~ 50/ias. Over a few hours, 
this corresponds to an expansion of ~ 1000km s _1 , consistent with the visibility 
modelling analysis. Most importantly, at 2 kpc an expansion in the plane-of- 
the-sky at lc would result in a closure phase rate of 180° in 8 minutes, which 
would be easily seen, especially in such high-SNR data. Though the phase 
rate is dependent on the degree of asymmetry of the structure change, there 
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Figure 3. Closure phase every 10 minutes from September 16 (dots) and 
17 (crosses) (Vivek Dhawan, priv. comm.) The closure phase from each day 
is offset by 20° for clarity . A number of different three antenna groups are 
shown, all with very high SNR data. 



is no evidence of relativistic expansion in the closure phase. This leads to the 
conclusion that the extensions are an artefact of calibration error. 



5. VSOP2 and stellar radio astrophysics 

Stellar observations with VSOP were particularly challenging due to the lack 
of the necessary sensitivity on HALCA to directly detect fringes. However, the 
observations of LS I +61° 303 provide some useful insight for planning stellar 
radio observations with future space- VLBI missions such as VSOP2. 

It is important that sensitivity is maximized. This implies the use of the 
most sensitive ground-based telescopes, not only in terms of collecting area (e.g. 
phased- VL A, Effelsberg, GBT etc), but also the observing frequency. This sug- 
gests 8.6 GHz will be the most suitable frequency since it offers an ideal combi- 
nation of sensitivity on the ground and, for non-thermal sources, a higher flux 
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than at the higher frequency observing bands specified for VSOP2. Ideally, the 
most sensitive 8.6-GHz system possible should be available on the spacecraft. 

An imperative for stellar observations in future missions is a source switch- 
ing or "phase-referencing" capability with the spacecraft. This ability was not 
available with HALCA. First, it is important to be able to involve the spacecraft 
in fringe-finding observations in order to establish the clock and orbit geometry 
offsets. If a source is sufficiently bright, phase-referencing observations can be 
done to utilise self-calibration to solve for time evolution of these parameters 
and the gain phase of the spacecraft. Further more, this provides an ability to 
improve orbit definition through astrometry. 

The VSOP observations of LS I +61° 303 were planned for a predicted 
time of outburst. Unfortunately, the onset of outburst was missed by 60 hours. 
Observations during outburst could only have been assured through a dynamic 
scheduling capability. Certainly if transient science is a key project for VSOP2 
then a dynamic scheduling capability that permits spacecraft command upload 
on short timescales is a necessity. 

Acknowledgments. I am very grateful to Vivek Dhawan (NRAO) for shar- 
ing his re-reduction of the ground-based d ata from the HALCA observation as 
part of a VLBA project on LS I +61° 303 ( Dhawan et al. 20061 ). Many thanks 
to Russ Taylor, Bill Scott a nd Marta Per e caula for reducing and analysing the 
data originally published in Taylor et al. ( 20001 ). And lastly, to Ed Fomalont 
and the VSOP team - many thanks for all the advice and help. 



References 

Albert, J., Aliu, E., Anderhub, H., et al. 2006, Science, 312, 1771 

Casares, J., Ribas, I., Paredes, J. M., et al. 2005, MNRAS, 360, 1105 

Dhawan, V., Mioduszewski, A., & Rupen, M. 2006, in VI Microquasar Workshop: Mi- 

croquasars and Beyond 
Fender, R. P., Stirling, A. M., Spencer, R. E., et al. 2006, MNRAS, 369, 603 
Frail, D. A. & Hjellming, R. M. 1991, AJ, 101, 2126 
Gregory, P. C. & Taylor, A. R. 1978, Nat, 272, 704 

Grundstrom, E. D., Caballero-Nieves, S. M., Gies, D. R., et al. 2007, ApJ, 656, 437 

Hermsen, W., Swanenburg, B. N, Bignami, G. F., et al. 1977, Nat, 269, 494 

Hutchings, J. B. & Crampton, D. 1981, PASP, 93, 486 

Johnston, S., Manchester, R. N., Lyne, A. G, et al. 1992, ApJ, 387, L37 

Maier, G. 2007, ArXiv e-prints, 709 

Massi, M., Ribo, M., Paredes, J. M., et al. 2004, A&A, 414, LI 

Massi, M., Ribo, M., Paredes, J. M., Peracaula, M., Estalella, R. 2001, A&A, 376, 217 
Paredes, J. M., Marti, J., Peracaula, M., & Ribo, M. 1997, A&A, 320, L25 
Paredes, J. M., Marziani, P., Marti, J., et al. 1994, A&A, 288, 519 
Peracaula, M., Gabuzda, D. G, & Taylor, A. R. 1998, A&A, 330, 612 
Ribo, M., Paredes, J. M., Moldon, J., Marti, J., & Massi, M. 2008, A&A, 481, 17 
Romero, G. E., Okazaki, A. T., Orellana, M., & Owocki, S. P. 2007, A&A, 474, 15 
Stirling, A. M., Spencer, R. E., de la Force, C. J., et al. 2001, MNRAS, 327, 1273 
Tavani, M. & Arons, J. 1997, ApJ, 477, 439 

Taylor, A. R., Dougherty, S. M., Scott, W. K., Peracaula, M., & Paredes, J. M. 2000, 
in Astrophysical Phenomena Revealed by Space VLBI, ed. H. Hirabayashi, P. G. 
Edwards, & D. W. Murphy, 223-226 

Taylor, A. R. & Gregory, P. C. 1984, ApJ, 283, 273 

Taylor, A. R., Kenny, H. T., Spencer, R. E., & Tzioumis, A. 1992, ApJ, 395, 268 



